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The ability to classify a state of matter from knowl- 
edge only of its many-body ground state is an exciting 
prospect in condensed matter physics. One established 
method is to evaluate the ground-state expectation val- 
ues A(j) = Gs\M(i)\^ Gs) for all local operators M(j). 
These order parameters distinguish the phase of matter 
in the Landau symmetry-breaking paradigm.^ Given a 
set of A(i), we ask: is the phase of matter uniquely de- 
termined? We now know it is not, an exceptional case 
being fractional quantum Hall states with, for example, 
filling factors z/ = 1/3 or = 1/5. These states are 
featureless insulators with all A(j) vanishing, yet they 
are distinguished by bulk topological properties which 
can be measured in Hall conductance experiments.'^ The 
discovery that symmetry-breaking does not classify all 
quantum phases has led to an explosion in the field of 
topological order.l^l A theory for the full classification of 
topologically ordered states is not yet understood. In 
this article we focus on one such classification method: 
the use of quantum entanglement to classify quantum 
ground states. 

Quantum entanglement is a recently fashionable ap- 
proach to classifying myriad states of matter ranging 
from spin-systems to topological insulators to fractional 
quantum Hall states.-^^^The primary classification tools 
have been the entanglement entropy, the topological 
entanglement entropy and the entanglement spectrum 
(ES); these tools are used to reveal subtle non-local 
correlations in topologically nontrivial systems. In this 
work, we specifically study the properties of the entangle- 
ment spectra of 2D time-reversal breakin^^sHUl ^nd time- 
reveral invarianlP2HlZ! topological insulators with spatial 
entanglement cuts. 

If the presence of edges does not break the symme- 
try that stabilizes the topological insulator {e.g. time- 
reversal, charge-conjugation, point-group symmetries), 
the physical spectrum must include gapless modes that 
separate the non-trivial insulator from the vacuum. Our 
formalism relies on topological insulators that manifest 
such edge states. There exists topological insulators 
which cannot be described by this formalism. A case 



in point is the topological insulator stabilized by inver- 
sion symmetry; since inversion symmetry is broken by 
the presence of edges, such insul ators d o not exhibit 
topologically-protected edge states.^ '^^ * ^'* * 

If these topologically-protected edge states are present, 
they interpolate across the energy gap and spectrally con- 
nect the conduction and valence bands - this is called 
spectral fiow. We will show that a subset of edge states 
at the Fermi level of the physical energy spectrum trans- 
lates to discontinuities in the trace of the single-particle 
ES, which we call a 'trace index'. We find that any free- 
fermion topological insulator that exhibits spectral flow 
has a non-vanishing trace index, which provides us with 
a new description of topological invariants. In this ar- 
ticle our case studies are the Chern insulator and the 
Z2 time-reversal invariant insulator. We investigate the 
trace indices in both case studies and demonstrate that 
the trace index is equivalent to the Chern number {Z2 
invariant) for the Chern {Z2) insulator. This formalism 
extends the work of Ref. |68l which presents a clear con- 
nection between spectral flow in the ES and the bulk 
topological invariant. 

In addition to the trace index, we identify more signa- 
tures of spectral flow in the single-particle and many- 
body ES; in this process we will establish a connec- 
tion between entanglement and quantum Hall physics. 
When an appropriate flux is threaded to induce charge 
transport, we flnd that the change in the full trace of 
the single-particle ES is exactly the quantity of charge 
that is pumped. Furthermore, we study insulators 
with global C/(l) symmetry (i.e. global charge conser- 
vation) so the many-body ES decouples into indepen- 
dent particle-number sectors, forming a 'particle-number 
space'. When charge is pumped in an adiabatic cycle, 
there is a global translation in the particle-number space 
by the quantity of pumped charge. In both case studies 
of the Chern and Z2 insulators, the amount of pumped 
charge is directly related to the topological invariant. 
Hence, these signatures in the single-particle and many- 
body ES separately provide us with new means to extract 
the Chern number and the Z2 invariant. 
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The outline of our paper is as follows. In Sec. [T] we 
introduce our chosen geometry, derive the single- and 
many-body entanglement spectra for free-fermion lattice 
models and introduce concepts such as an entanglement 
ground state and an entanglement Fermi level. Our case 
studies are the Chern insulator and the Z2 time-reversal 
invariant insulator (also called the quantum spin Hall 
insulator). In Section |ll] we define the U(\) trace in- 
dex for the Chern insulator and show its equivalence to 
the Chern number for translationally invariant and dis- 
ordered systems. In Section |III| we define the trace 
index for the insulator and show its equivalence to the 
Z2 invariant. In addition, we identify the signatures of 
spectral flow in the single- and many-body entanglement 
spectra of both case studies and establish a connection 
between quantum Hall physics and entanglement. We 
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make some concluding remarks in Section IV 



I. ENTANGLEMENT SPECTRUM IN 
FREE-FERMION LATTICE MODELS 

Given a many-body ground state I^'gs) in a Hilbert 
space C and its corresponding density operator p = 
I^Gs) (^Gsl: we partition the Hilbert space into two parts 
C — A® B. The reduced density matrix in the sub- 
system A may be expressed as the exponential of an en- 
tanglement Hamiltonian, A/'e"^""', which is normalized 
such that Tr/3_4 = 1. This is likened to the density oper- 
ator of a thermodynamic system with Hamiltonian 
and unit temperature.!^ 

In this article we distinguish between the physical 
ground state I^'gs) and the entanglement ground state 
15*^3), which is the lowest-energy many-body eigenstate 
of the entanglement Hamiltonian H^^-t. The limit of zero 
entanglement is realized if the many-body ground state 
|4'gs) is a single product state in the two subsystems. In 
this case the entanglement energy gap between l^'^s) and 
its excited states is infinitely large and the entanglement 
entropy is zero. A highly entangled I^'gs) is character- 
ized by small entanglement-energy excitations of |^'g"|)- 

For a topological insulator, the entanglement entropy 
of at least one partition A^B cannot be tuned to zero by 
any change of the physical Hamitonian H that preserves 
the gap in the physical spectrum and the symmetry that 
stabilizes the topological phase. This implies that the en- 
tanglement energy gap between the entanglement ground 
state and its excited states cannot be tuned to infinity.!^ 
The topology of the phase is thus encoded in the spec- 
trum of Pa, which we call the many-body entanglement 
spectrum (ES). 

For free-fermion Hamiltonians, the many-body ES may 
be derived from t he eig envalues of a reduced one-body 
correlation matrijP^I^, which we henceforth call the 
single-particle ES. The reduced correlation matrix is de- 
fined in Section [l A| and its connection to the many-body 
ES is derived in Section [TB} In addition, we discuss the 
differences between the many-body entanglement spectra 




FIG. 1. (a) Our chosen geometry is a square lattice with open 
boundary conditions along x and periodic boundary condi- 
tions along y: an open-end hollow cylinder. Quantum Hall 
phases such as the Chern insulator exhibit modes that extend 
around the edges of the cylinder. Extended states on opposite 
edges have opposite chirality; the direction of propagation of 
these modes are indicated by red and blue arrows, (b) We 
perform a spatial entanglement cut by partitioning the cylin- 
der into two halves - the left half is labelled A and the right B. 
For the Chern insulator, the spectrum of the reduced correla- 
tion matrix in region A includes edge modes (green arrow) 
on the entanglement cut between A and B. These edge modes 
contribute maximally to the entanglement entropy. 



of a trivial and topological insulator in Section |IB 



Single-particle entanglement spectrum in 
free-fermion lattice models 



In this section we review the construction of the single- 
particle entanglement spectrum (ES) from a free-fermion 
lattice Hamiltonian H. We focus on 2D lattice models: 



H 



'Hj Ha'-jP 



(1) 



with spatial indices i,j and orbital/spin indices a, {3. 

We will be discussing topological insulators which 
manifest edge states when placed on a lattice with open 
boundary conditions. We work with a cylindrical geom- 
etry so that the physical spectrum includes states that 
extend around the edges of the cylinder. We choose ij 
to be the azimuthal direction and x to lie parallel to the 
symmetry axis. We have open boundary conditions at 
sites X = 1 and x = N^, as shown in Fig. [T]-a. 

If the Hamiltonian is translationally- invariant, we 
Fourier transform equation ([T]) in y and obtain 



H 



E.4. 

ky ^i. j 



j0 



(2) 



with i,j now labelling the x-coordinates of lattice sites. 
We sum over repeated indices and express the n-th single- 
particle eigenstate of (|2j) in Bloch form: {y\^ky) = 
g*fc!/y[yn^]a Ijq,) with \ja) a basis state with orbital/spin 
a and i-coordinate j. 

The Hamiltonian decouples into irreducible represen- 
tations labelled by the conserved momentum ky] in 
each representation, the n-th eigenstate of the Bloch 
Hamiltonian hf^^ {ky) has the corresponding projector 
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inrfiky) = KJfK:t The many-body ground state 
wavefunction I^'gs) is a single Slater determinant of all 
single-particle eigenstates of H with energies less than the 
Fermi energy \^as) = n„,fc„;e„(fc„)<,. 7^ JO) - ^ith 
normal mode operators -/Iky = KJj/sclyji^- 

We define a one-body correlation matrix G which is ex- 
pressible as the complex conjugate of a sum of projectors 
of occupied single-particle states: 



[G]tfiky) 



kyia'^kyjP I^C 



(3) 



n;en(ky)<^i 



Being a sum of projectors, the correlation matrix must it- 
self be a projector, i.e. G{ky)^ = G{ky), with eigenvalues 
or 1. 

To make a spatial entanglement cut, we partition the 
cylinder into two halves A and B each of length Nx/1\ 
this is illustrated in Fig. [T]-b. We choose the conven- 
tion that region A is the left half of the cylinder. We 
define the reduced one-body correlation matrix G'^{ky) 
as the correlation matrix G{ky) with its spatial indices 
restricted to the subset of sites in region A, 



'^kyia^kyj ^ I 



(4) 



If there are Ni^^^i local degrees of freedom per site, G^{ky) 
has = NxNi„^^i/2 eigenvalues which we label by 
Ci{ky). The spectrum of G'^lky) is bounded by and 
1. If the Hamiltonian is completely on-site, i.e. diago- 
nal in real space, the spectrum consists only of O's and 
I's and consequently the entanglement entropy is zero. 
This is the atomic limit of an insulator. We are inter- 
ested in topological insulators which are defined to be 
not adiabatically connected to the atomic limit; these 
phases have finite entanglement entropy which cannot be 
removed unless the insulator gap collapses or the symme- 
try that stabilizes the insulator is removed. 

Eventually we will describe disordered systems, so we 
want to define the single-particle ES when there is no 
translational symmetry in y. In this case, the correla- 
tion matrices G and G^ do not decouple into different ky 
sectors so the generalization is implemented by removing 
the ky labels in the definitions in Eq. ([3| and 



B. Relationship between the many-body and 
single-particle entanglement spectra of free-fermion 
models 

In this section we derive a relation between the many- 
body entanglement spectrum (ES) and the spectrum of 
the reduced correlation matrix G^ and introduce the con- 
cepts of an entanglement ground state and an entangle- 
ment Fermi energy. We illustrate these ideas with two 
examples - the entanglement spectra of a trivial and topo- 
logical insulator - and discuss their differences. 



The many-body ES is the set of eigenvalues of the the 
reduced density operator in region A. If we have trans- 
lational symmetry in y, we may factorize Pa into distinct 
ky sectors and express it as the exponential of an entan- 
glement Hamiltonian iJe„t: 



Pa 



^PA{ky)^N\{ 



^-H^nt{ky)i 



(5) 



In the rest of the article we let the indices i,j of the 
matrices G,G^, iJ„„t denote collectively all the degrees of 
freedom: spatial, spin and/or orbital. 

The matrix H^^^(ky) has Qa eigenvalues - the single- 
particle entanglement energies - which we label by S,i{ky); 
ni{ky) are defined as the occupation numbers of the cor- 
responding single-particle eigenstates. Af normalizes the 
reduced density operator such that it has unit trace; in 
the thermodynamics analogy J\f~^ is the partition func- 
tion of an entanglement Hamiltonian with unit tempera- 
ture. In this analogy, the final form of the reduced density 
operator assumes a familiar form 



PA{ky) = Y[ 



i{ky)X 



1 + e^^'C^B^ 



(6) 



with normal mode operators xl^^i = Vj{ky)cl^j. 

Let us define Tr . ^ as the trace over all states in sub- 
system A in the momentum sector ky. We note that 
TrA,fcyPA(fc«) = 1- The reduced correlation matrix as 
defined in El) can be expressed with reduced density op- 



„t 



erators as Gf^(fc„) = Tr^,.^ 
Appendix [K\ the final steps to derive 



[G^iky 



1 

1 -|- gffont(fci,) ■ 



We show in 



(7) 



which relates the single-particle entanglement energies 
(^i) to the eigenvalues (Ci) of G^. Given a set of oc- 
cupation numbers {ni{ky)} for the single-particle en- 
tanglement eigenstates, the corresponding eigenvalue 
A[{ni(fcj,)}] of the reduced density operator is 



A {n:{ky)} ^liUil-Qiky) 

' ky i=l ^ 

= n n c^iky) n 

ky 2 G occupied J £ unoccupied 



ni(ky) 



(8) 



With this relation we can construct the many-body ES 
(i.e. all eigenvalues of Pa) by occupying the single- 
particle entanglement states for each ky in all possible 
combinations. 

The task of finding structure in the many-body ES is 
eased by exploiting the global U{1) charge conservation 
symmetry of the Hamiltonian to decouple the many-body 
reduced density matrix into different particle number sec- 
tors. Let us clarify what this means. The density matrix 
over the entire Hilbert space is the projector \^gs) (^'gsI 
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FIG. 2. Top: Figures |2]-a and b respectively show the entan- 
glement energy spectrum and the many-body entanglement 
spectrum (ES) of a trivial insulator that is adiabatically con- 
nected to the atomic insulator. The latter is scaled logarith- 
mically. The many-body ES belongs in the half-filled particle 
number sector, i.e. half the total number of single-particle 
entanglement eigenstates are occupied. Bottom: Figures [2]-c 
and d respectively show the same spectra for a topological 
insulator which exhibits spectral flow. The many-body ES is 
plotted as a function of the total momenta ky of the many- 
body entanglement state. As a reference point, we define ky^ 
as the total momentum of the entanglement ground state. 



of the ground-state many-body wavefunction, which has 
one particle number Nqs - the number of occupied single- 
particle states in the ground-state Slater determinant. 
When we partition the cylinder we have a finite proba- 
bility of trapping any number of particles in region A; we 
are only limited by the original number of particles we 
began with, and how many fermionic particles can fit in 
A. We are often interested in many-body entanglement 
eigenvalues in a single particle number sector - N^. In 
this case, we fill up the single-particle entanglement levels 
using every possible combination of occupation numbers 
{ni{ky)} that sum up to N^^ - this is the set of many-body 
entanglement eigenvalues in the sector. We define the 
set of all nonzero eigenvalues in the Nj^ sector as 

{4 ={n n c.(fc.) n 

Na ^ ky i^occupied j^unoccupicd 

J2'l.{ky) ^ (9) 

In analogy with the physical Hamiltonian, we may con- 
struct the entanglement ground state in the sector 
by filling up of the lowest-energy eigenstates of the 
entanglement Hamiltonian. This ground state has a one- 
to-one correspondence with the largest many-body eigen- 



value A in the sector through relation ([s]). The en- 
tanglement Fermi energy is the largest single-particle 
entanglement energy in the many-body ground state. All 
other eigenvalues A in the sector correspond to ex- 
cited many-body states which derive from particle-hole 
excitations of the entanglement ground state across this 
entanglement Fermi level. 

As an example, we plot in Fig. [2]-a the single- 
particle entanglement energy spectrum {^} of a trivial 
translationally-invariant insulator that is adiabatically 
connected to the atomic insulator. There is a clear spec- 
tral gap which separates the lowest positive-^ band and 
the highest negative-^ band. We consider the half-filled 
entanglement ground state with all the negative-energy 
single-particle entanglement states filled. We denote the 
total momenta of all single-particle states in this ground 
state as k^^ . In Fig. [2|b we plot the many-body ES {A} 
of the half-filled sector as a function of total momentum 
ky] more precisely, we have plotted {—In A} with A the 
eigenvalues of the reduced density matrix in region A. 
Through relations ([t]) and (|8]), the entanglement ground 
state is in one-to-one correspondence with the smallest 
value of —In A in the many-body spectrum at momen- 
tum ky — ky^ . All other values of A in the half-filled 
sector correspond to excited many-body entanglement 
states which derive from particle-hole excitations across 
the entanglement gap. Since these particle hole excita- 
tions are gapped in the single-particle spectrum of {<^}, 
there is also a clear spectral gap in {—In A} separating 
the the many-body ground state from its excited states. 

We compare these plots with that of a topological in- 
sulator in Fig. |2]-c and d. In the next Section we iden- 
tify these new plots with a Chern insulator with Chern 
number —1. As shown in Fig. [2]-c, there is a set of 
single-particle entanglement eigenstates that interpolate 
across the gap and connects bands with negative and 
positive ^. This property is called spectral flow. These 
single-particle states in the gap extend over the spatial 
entanglement cut between subsystems A and B, so we re- 
fer to them collectively as the entanglement edge mode. 
The entanglement ground state in the half-filled sector 
has a Fermi level S^p — Q which intersects the dispersion 
of the entanglement edge states. Thus, the low-energy 
particle-hole excitations of the ground state are gapless 
and correspond to excitations of the entanglement edge 
mode. This is reflected in the many-body spectra (Fig. 
[2]-d) by the absence of a spectral gap between the entan- 
glement ground state and its excited states. In the next 
Section we present our flrst case study of a topological 
insulator with spectral flow - the Chern insulator. 



II. CHERN INSULATOR 

The Chern insulator (CI) realizes a quantum anoma- 
lous Hall effect without application of a net magnetic 
fleld.*® The simplest CI model exists on a square lattice 
with two degrees of freedom per site and has the Hamil- 
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FIG. 3. Top: Figures 2-a,b,c show, respectively, the physical energy spectrum, the single-particle entanglement spectrum (ES), 
and the trace of G^{ky) for a trivial insulator (Ci = 0) on an open-end cylinder. In Fig. 2-a and b, bulk bands on opposite 
sides of the gap are spectrally disconnected from each other. In Fig. 2-c, the trace of is a continuous function of ky so the 
trace index Au{i) ~ 0. Bottom: Figures 2-d,e,f show, respectively, the physical energy spectrum, the single-particle ES, and 
the trace of G^{ky) for a nontrivial insulator with Ci = —1. In Fig. 2-d and e, a set of edge states interpolate between bulk 
bands on opposite sides of the gap - this spectral flow distinguishes trivial and nontrivial Chern insulators. In Fig. 2-f, the 
trace of G^ drops discontinuously by 1 at fc^ = - the trace index Au{i} ~ ^1- 



tonian 

H = 



E 

t' 



(4+1,1 



+ (2 - m) {clyO^c^^y)] , (10) 

where Cx,y = [cx.yf (^x,yi) destroys a spin up/down 
fermion on site (x, y), and t,t' ,m are model parameters. 
For simplicity we set t = t' = 1 and the phases of the 
model are controlled by m alone. Full translational sym- 
metry allows us to express the Hamiltonian as 

H = ^^c^ [sin kx (Ji + sin ky ai -f M{k) CT3] c^. 



M{k) = 2 — m — cos — cos ky 



(11) 



For m < and to > 4 the insulator is trivial and has 
Chern number Ci = 0. For < to < 2 the insulator has 
Ci = — 1 and for 2 < TO < 4 the insulator has Ci = +1. 
To illustrate the different nature of these phases, we plot 
the physical energy spectra for the Ci = and Ci = — 1 
cases with our chosen geometry in Figs. |3[a and d re- 
spectively. The two spectra are distinguished by gapless 



chiral edge states which produce a quantum Hall effect 
in the nontrivial phase. 



Let us now examine the single-particle entanglement 
spectrum of the CI. The eigenvalues of the full corre- 
lation matrix G are either or 1 since G is a projector. 
With an entanglement cut, most of the eigenvalues of the 
reduced correlation matrix G"* lie close to or 1. For the 
trivial insulator with Gi =0, the bands of eigenvalues 
near 1 are spectrally disconnected from the bands near 

0, as shown in Fig. [3]-b. If the insulator is nontrivial, 
there is in addition a set of interpolating single-particle 
entanglement states that flow across the gap from to 

1, as shown in Fig. [3]-e. This set of states extend over 
the entanglement cut between regions A and B, and con- 
sequently contribute the most to entanglement entropy. 
The presence of spectral flow is a topological property 
that guarantees that the entanglement entropy is nonzero 
for spatial cuts. 
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Single-particle Entanglement Spectrum 
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FIG. 4. Left side: Figures 3-a,c,e,g show schematically the physical energy spectrum, the single-particle entanglement spectrum 
(ES), the trace of the single-particle ES and the entanglement energy spectrum respectively, as functions of ky, for a many-body 
ground state of a nontrivial Chern insulator. Right-side: Figures 3-b,d,f,h label the new functions when a flux quantum is 
inserted through the symmetry axis of the cylinder. These are schematic flgures that emphasize importance features in the 
physical and entanglement energy gaps; bulk bands away from the gaps are conveniently portrayed as dispersionless. Eigenvalues 
that are colored red, green and blue correspond to eigenstates that extend over the left edge, the entanglement cut and the 
right edge of the cylinder respectively. Occupied single-particle states of a many-body Slater determinant are labelled by black 
circles in both the physical and entanglement energy spectra. 



A. CI: The trace index and its equivalence to the first discuss a clean, free-fcrmion insulator. We illustrate 
Chern number the index for a CI with Ci = -1, Fermi energy ^ = and 

edge states crossing at ky = as in Fig. |4]-a. We consider 

In this section we define the trace index for a CI and 
demonstrate its equivalence to the Chern number. We 
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the trace of the reduced correlation matrix G^{ky): 



Tr [G-{ky)\ = (vf, 



IE 4- 



|*GS) . (12) 



as the number op- 



We identify NA{ky) = E 
erator in region A and in the ky sector, so we estabhsh 
that the trace of G^{ky) is the ground-state expectration 
value of NA{ky). This suggests that the trace is a contin- 
uous function of ky except where modes cross the Fermi 
level. Since we have assumed we are in a bulk gap, modes 
crossing the Fermi-level have to be edge states. In this 
example, one such mode, colored red in Fig. |4]-a, crosses 
the Fermi level between ky = and jf- . To the immediate 
left of the crossing at ky = 0, the occupied state just un- 
der the Fermi level is on the left edge (up to exponential 
finite-size corrections) and contributes -1-1 to the trace of 
0^(0); to the immediate right at ky ^ j^, the occupied 
state just under the Fermi level is on the right edge and 

Hence there 



does not contribute to the trace of G'*(^). 



-1 which we 



is a discontinuity: TrG^(|^) - TrG^(O) 
tentatively identify as the Chern number. 

We now consider the most general insulator. The 
Fermi energy may take on any value in the gap that does 
not intersect an extended state; this keeps the phase in- 
sulating. We define i^orosa as the average of two adjacent 
discrete momenta where a left-edge dispersion intersects 



the Fermi level; in the above example, -ftTcroas — 



If 



the gradient of the edge dispersion at such an intersec- 
tion is positive {i.e. ^^(i^croBa) > 0), there is at least 
one occupied left-edge state to the left of the crossing at 
^croBB — We claim that this occupied left-edge state 
is an eigenstate of G'^iK^,^^^ — with eigenvalue ex- 
ponentially close to 1, so it contributes 1 to the trace of 
G^. We claim that the unoccupied left-edge state just 
above the Fermi level is an eigenstate of G"^ (K^,^^^ + ) 
with eigenvalue exponentially close to 0, so it does not 
contribute to the trace of G^. As a consequence, the 
trace of G"* decreases discontinuously by 1 across this 
crossing. An analogous argument leads to the conclu- 
sion that a crossing with a negative gradient results in a 
discontinuous increase by 1. 

These claims are justified by the following theorem. 

Theorem II. 1 As defined in equation Q), let G be the 
one-body correlation matrix of a many-body eigenstate of 
a free fermionic Hamiltonian H . Let the Hilbert space 
be partitioned into a direct product ^ (g) _B. Let m^^^ be a 
single-particle eigenstate of the Hamiltonian that is con- 
fined in A. As defined in equation Q), let G* be the one- 
body reduced correlation matrix with indices restricted to 
A. If the many-body wavefunction includes the single- 
particle state u^-^\ then u^^-** is an eigenstate of G^ with 
eigenvalue 1. If the many-body wavefunction excludes the 
single-particle state u^^\ then u'--^^* is an eigenstate of 
G^ with eigenvalue 0. 



The theorem has the following physical motivation: a 
state that is localized on the physical edge of A has ex- 
ponentially vanishing probability to be in region B. In 
the thermodynamic limit, this edge state is fully confined 
in A and cannot entangle regions A and B. As per our 
discussion in Sec. |IB[ the single-particle states that do 
not contribute to entanglement entropy have eigenvalues 
or 1 in the spectrum of G"*. The proof of this theorem 
is presented in Appendix [Til 

We apply this theorem to our example of Fig. |4]-a, 
where we show that the occupied left-edge state at fc^ = 
is an eigenstate of G^{Q) with eigenvalue 1, also colored 
red in Fig. jijc. The vacant left-edge state At ky — jf- 
in Fig. |4]-a is shown to be a zero-eigenvalue eigenstate of 
G'^(;^) in Fig. ijc. In Fig. |4|e we plot the trace of G^ 
and show that the discontinuity in the trace lies between 
ky — and where the left-edge mode intersects the 
Fermi level. 

More generally, there may be n(+) (n(_)) left-edge 
states crossing the Fermi energy at a given crossing K„„^^ 
with positive (negative) and each state contributes 
— 1(-|-1) to the trace discontinuity. The net discontinuity 
is «(„) — We define the U{1) trace index {Au(i)) 

for clean, free-fermion insulators as the sum of the trace 
discontinuties over all K„g^^ in the thermodynamic limit: 



-^-.-.3- ETr[G^](i^... + ^) 



^Tr[Gl(i._--) 
= "(-)(^=.o.)-n,+,(if_) (13) 

The index effectively counts the net number of chiral 
single-particle states, on one edge of the cylinder, that 
intersect the Fermi level. The edge states are responsi- 
ble for the quantum anomalous Hall effect when fiux is 
inserted through the symmetry axis of the cylinder, or 
equivalently when an azimuthal electric field is applied. 
•^u(i) is precisely the number of charges pumped from one 
edge to another during an adiabatic insertion of a unit 
fiux, i.e. ji = e^^ j^EjAuii)- This connection between the 
Hall current and the trace index is elaborated in Section 
In summary, we have identified the U{1) trace 



IIBl 



index with the first Chern number (Gi) of the CI. 

-4(7(1) — Cl. 



(14) 



As a confirmation of ( 14 ) , we numerically plot the trace 
of G-* for a trivial (Gi — 0) and nontrivial (Gi = — 1) 
insulator respectively in Fig. [Sj-c and f. The trace is a 
continuous function of ky for the trivial case, and drops 
discontinuously by 1 at fcj, = for the nontrivial case. 

We remark that Au(i) is quantized to an integer in 
the thermodynamic limit. The first reason is that physi- 
cal edge states are perfectly confined in region A in this 
limit and this is a necessary condition to apply Theo- 
rem II. 1 for any finite-size lattice the trace index has 
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exponentially small deviations from an integer because 
an edge state predominantly in A has exponentially van- 
ishing probability to be found in B. The second reason is 
that physical bulk states do not contribute to the index 
in this limit. While states in the valence band have finite 
probability to be in A and thus contribute to the trace of 
G^, the projectors of the bulk states are a smooth func- 
tion of ky in the limit that the discrete spacing goes 
to zero. Any discontinuity must arise from a difference 
in occupation numbers between two adjacent bulk states 
- this is impossible if the Fermi energy is in the bulk gap. 



B. CI: Flux threading and entanglement 

For a nontrivial CI, a time-dependent flux induces an 
azimuthal electric field which drives a quantum Hall cur- 
rent. We investigate the effect of flux threading on the 
single-particle entanglement spectrum in Section |II B 1| 
and establish a connection between entanglement and 
quantum Hall physics. In Section |IIB 2| we identify the 
signature of spectral flow in the many-body entanglement 
spectrum as a global translation by Ci of all particle- 
number sectors. In Section [ll B 3| we generalize the trace 
index to describe disordered CI's by replacing ky with 
flux as a pseudo-momenturrPSl. In addition we show that 
a disorder-induced topological phase transition occurs 
when the disorder strength is increased beyond a criti- 
cal value. 



1. Effect on the single-particle entanglement spectrum 



Adiabatically inserting a unit flux $o = f into a non- 
trivial CI leads to an excited many-body state. The 
flux insertion is implemented by a Peierls's phase which 
may be gauged away at the cost of twisting the bound- 
ary conditions. This twist changes the allowed values 
of momenta ky - each single particle state transforms as 

fcy + F^^^^ where $(f) is the instantaneous 
flux at time t. Upon insertion of one flux quantum, the 
physical energy spectrum is invariant but there is a shuf- 
fling of occupation numbers as single-particle states at ky 
transform into adjacent states at fcj^ -I- 

We illustrate this process in Fig. |4j Before flux- 
insertion, the Fermi energy lies over a crossing of left- 



\ky) 



and right-edge mode dispersions at K^^^ 



as shown 



in Fig. |4]-a. After flux-insertion, the resultant Slater de- 
terminant is shown in Fig. |4]-b. Since the hybridization 
between states on opposite edges is exponentially small, 
the left-edge state at ky = transforms, upon flux inser- 
tion, into a left-edge state at ^ . The right-edge state at 

ky = ^ transforms into a right-edge state at ky = 
These transformations lead to a particle-hole excitation 
of the ground-state and the resultant wavefunction takes 



operators of the lowest energy on the left and right edges 
respectively. A net charge equal to the Chern number 
(—1 in this example) is pumped from the left to the right 
edge; this is a consequence of spectral flow in the physical 
energy spectrum. 

Let G;* and G-i denote the reduced correlation matri- 
ces before and after insertion respectively; their spectra 
are plotted in Figs. |4]-c and d respectively. The dif- 
ference between the two matrices is that G^ has one 
extra unit eigenvalue and one less zero eigenvalue at 
ky — j^. This difference arises because the flux-inserted 

state has an extra occupied left-edge state at fcj^ = 
We observe in Fig. [4[e that the ground-state satisfies 



Tr G;^ (1^ ) - Tr G;;^ (0) = - 1 , whereas in the flux-evolved 



state, TrG^^(^) - TrG^^_^(^) = -1 as shown in Fig. 
|4|f. The position of the trace discontinuity at K^^^^^ = 
is translated by ^ when a unit flux is inserted, but the 
magnitude of the trace discontinuity and the U{1) trace 
index Au(i) is invariant. 

More generally, there may be multiple crossings K„^^^ 
of the left-edge states with the Fermi level, and each 
crossing may contribute ±1 to the trace index depend- 
ing on the gradients -^{K^^^^^). When an integer flux is 
inserted, all the positions of the discontinuities are trans- 
lated by i^cross — > J'^cross + If" tut Au(i) Is iuvariaut. We 
leave the proof of this claim to Appendix |III| In Section 
we identified the A,,n-\ with the Chern number of the 



ground state; in this Section we show that this equality 
may be extended to any flux-inserted state which only 
differs from the ground state in the occupation numbers 
of its edge states. This implies that we may separately 
extract the Chern number from both the ground-state 
and its flux- inserted state through Au(i)- 

We now present an equivalent way of extracting the 
Chern number from two many-body states differing 
in fiux. While the physics is identical, we focus on 
a different quantity in the entanglement spectrum to 
emphasize the connection with the quantum Hall ef- 
fect. This quantity is the full trace of G^ over all ky 
sectors Tr G-{ky) = (^| J*) = 

{'i'\NA 1^)- This is the expectation value of the total 
number of particles in region A. We claim that 

^Tr K(fc,)] -5]Tr [G:iky)] = -C,. (15) 



the form |5'$„ 



Proof: under flux insertion, the positions of the discon- 
tinuities in the trace of G'^[ky) are translated by ^ but 
the magnitudes of the discontinuities are invariant. This 
implies that each discontinuity of — 1(-|-1) when trans- 
lated forward adds (deducts) ^ to the area under the 
piecewise-continuous function G'*(fcy). The net decrease 
in the area under G^{ky) is the sum of the discontinu- 
ities. The full trace of G^ is just the area under G^{ky) in 
units of "j^ . By definition, the sum of the discontinuities 



^'gs) with ^i^,"fji normal mode is the index Aun) which is also the Chern number. 
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We remark that Tr [G^] and the charge pumped from 
B to A during an adiabatic cycle are well-defined quan- 
tities independent of the details of the Hamiltonian. In 
disordered and/or interacting systems, ky is not a good 
quantum number but Tr [G^J - Tr \G^] = -Ci still 
applies. We substantiate by reproducing ( 15 ) through a 



quantum Hall response. The time-dependent flux gener- 
ates an azimuthal electric field 



En 



dAy 

dt 



1 a$ 



(16) 



For any CI (interacting/free, disordered/clean), a quan- 
tum Hall current flows in response to the field Ey: 



d (^\ N, I*) 



d Tr [G^ 



dt 



N„ 



dt 



(17) 



This current is shown to be proportional to the rate 
of change of the full trace of G^. This implies that 
the change in flux $ is proportional to the change in 
Tr [G'^] . When we insert one flux quantum ^, the change 
in Tr [G-*] equals -^cr„ = -Gi. 



2. Spectral flow and the many-body entanglement spectrum 

How spectral flow is manifested in the physical energy 
spectrum is well-understood - the global spectrum is in- 
variant under insertion of a unit flux but Gi charge is 
pumped resulting in a quantum Hall current. In this 
section we investigate how spectral flow manifests it- 
self in the many-body entanglement spectrum (ES). We 
study insulators with global U{1) charge conservation 
symmetry so the many-body ES decouples into indepen- 
dent particle-number sectors, forming a 'particle-number 
space'. We find that, analogously, the many-body ES is 
invariant under insertion of a unit flux and the charge 
pumping is manifested in a global translation of Gi units 
in particle-number space. 

We choose the convention that Gi charge is pumped 
from region A to B when a unit flux is inserted, so there 
are Gi less occupied left-edge states in region A. Theo- 
rem [TlT] informs us that in the thermodynamic limit the 
single-particle ES now has Gi less unit eigenvalues and 
Gi more zero eigenvalues. Let us illustrate this with the 
same example of a Gi = — 1 insulator with a Fermi en- 
ergy (/J, = 0) that intersects the left-edge mode dispersion 
at K^^„^^ — For this example, the spectra of the en- 
tanglement Hamiltonians before and after flux insertion 
are plotted in Figs. |4]-g and h respectively. The spectra 
in Figs. |4]-g and h are respectively related to the spectra 
of G"* in Figs, [ij-c and d through the relation ([7|. We 
first consider the A^^-number sector of the many-body ES 
before flux insertion. As illustrated in Fig. |4]-g, the en- 
tanglement ground state of this sector is constructed by 
filling up of the lowest-energy eigenstates of the en- 
tanglement Hamiltonian; we call the largest energy of the 
occupied states the entanglement Fermi energy ^p. The 



many-body entanglement spectra in this number sector is 
derived by applying Eq. ^ to the entanglement ground- 
state and all possible particle-hole excitations across the 
entanglement Fermi level. We denote the set of nonzero 
eigenvalues in this sector by {A(<I> = 0)}Ar^;^p. 

When a unit flux is threaded, a right-edge state is 
transported to the left edge. The difference between Fig. 
(before) and d (after) is that a eigenvalue of G^ at 
is raised across the entanglement gap and now has 

value 1. Equation (m) informs us that C(;^) — 1 corre- 
sponds to an entanglement energy ^(|f^) = — oo. In Fig. 

h (after) we show that the entanglement ground state 
formed from occupying all single-particle states under the 
same entanglement Fermi level now has an additional 
particle with ^ = — oo at fcj, = 

We ask what is the spectrum obtained from particle- 
hole excitations across the same Fermi level but now 
with Nj^ + 1 particles? In other words, what is {A(<I> = 
$o)}AfA+i;CF? We claim that it is identical with {A(<I> = 
0)}7Va:Cf the thermodynamic limit. 

The physical argument is as follows: we have shown 
that single-particle eigenstates of the physical Hamilto- 
nian that are confined in A cannot entangle regions A 
and B, thus they have eigenvalues = or 1. Hav- 
ing more or less of these zero-entanglement states as 
a result of flux insertion cannot affect the entangle- 
ment entropy. Since the two sets of nonzero eigenvalues, 
{A($ = 0)}Ar^;^p and {A($ = $o)}wa+1;Cf> differ only 
by one zero-entanglement single-particle state at the 
two sets must be identical for the entanglement entropy 
to be invariant. Let us define the many-body wavefunc- 
tions before and after insertion of the unit flux as |5'o) 
and respectively. The number sectors before and 

after insertion are related by 



(18) 



A ^ A 

The mathematical proof is presented in Appendix |IV| 
This is a global translation by Gi in particle-number 
space, i.e. all particle-number sectors flow to Na — Ci, 
which is analogous to the charge pumping in the physical 
energy spectrum. This is a change in the labels of each 
sector; the eigenvalues in each sector are invariant. As 
a consequence, the entanglement entropy, which depends 
only on the values of Ai, is invariant under integer flux 
insertion. 



The derivation of (18) through Theorem II. 1 and the 
proof in Appendix |IV| is applicable to any free-fermion 
Hamiltonian, so the particle-number flow occurs for dis- 
ordered Chern insulators as well. If we introduce inter- 
actions, this derivation no longer applies but the prin- 
ciple remains: with regard to the many-body entangle- 
ment eigenvalues, the many-body ground state with Gi 
less/more particles on the edge of A is indistinguishable 
from the many-body ground state before flux insertion, 
because the difference in the two wavefunctions lies far 
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away from the entanglement cut between A and B. If 
the interactions preserve the global U{1) symmetry, the 
many-body ES still decouples into particle-number sec- 
tors which flow as in ( 18 1. Wc also remark that equation 
( 18 1 is derived for the integer quantum Hall effect. How- 



ever, spectral flow is generic in all quantum Hall states 
and we have an analogous result for the fractional quan- 
tum Hall states if we insert more than one unit flux in 
analogy with the charge pumping picture presented in 
Ref. [m 



3. The trace index of disordered systems and a 
disorder-induced topological phase transition 

The effect of disord er on topological insulators is an on- 
going field of resarch.ti^l^2HZZI Generic topological insula- 
tors are characterized by extended bulk states that resist 
Anderson localization; these insulators are robust against 
disorder provided some symmetries are preserved.'^ For 
the nontrivial CI on a cylinder, we introduce disorder by 
adding a random spin-dependent potential on each site. 
The disordered Hamiltonian takes the form 



W 



x,y,a 



(19) 



with Ho the lattice Hamiltonian introduced in Eq. ( 10 1 



a denotes a spin index and Wx^y^a is a random number 
between —0.5 and -1-0.5. We tune the strength of this 
disorder potential by a parameter W. 



The trace index Au(i) as defined in (13) for clean, free 
systems relies on the conserved quantity ky. However, 



we have shown in ( 13 ) that Au{i) effectively counts the 
net number of chiral single-particle states, extended over 
either edge of the cylinder, that intersect the Fermi level. 
To generalize the description of trace indices to disor- 
dered free-fermion topological insulators, we replace ky 
with an adiabatic parameter - the flux inserted through 
the symmetry axis of the cylinder - that acts as a pseudo- 
momentum. 

The following discussion mirrors the introduction of 
Au{i) in Section II A and is conceptually similar. In anal- 



ogy with the Fermi energy of a clean insulator, we choose 
an energy e J in the mobility gap of the disordered energy 
spectrum. For each value of the flux $, we deflne the 
many-body wavefunction |$) as the Slater determinant 
of all single-particle states with energies less than £*: 
I*) = n„;e„($)<£| |0)- If the insulator exhibits spec- 
tral flow, inserting integer flux into some initial ground 
state 1$ = 0) creates particle-hole excitations. In this 
case, 1$), which is defined as a fully-filled Fermi sea for 
all values of the flux, is not necessarily the flux-evolved 
state that is adiabatically connected to |4> = 0), i.e. we 
are specifying for |<I>) that all of the lowest-energy single- 
particle states are filled regardless of how they adiabati- 
cally evolve. The reduced correlation matrix defined as 



[Cl"^($)= ($1 



1$); iJeA 



(20) 



is analogous to but different from our previous definition 
of in (|4|. i,j are lattice sites with the x-coordinate 
restricted to he in the range [l,Nx/2]. 

We define the C/(l) trace index (y^*^)) for disordered 
insulators as the sum of the discontinuities in the trace 
of as a function of the fiux $ in the range [0,27r). 
These discontinuities arise when states in region A dis- 
perse with $ and intersect e% . There exist localized states 
in the mobility gap which are confined in A but do not 
extend around the circumference of the cylinder. The 
effect of flux on such a state is a change of gauge: its 
wavefunction acquires a phase but its energy is invari- 
ant. These localized states do not disperse as a function 
of i> and cannot intersect e%. If the insulator is topologi- 
cal and exhibits spectral fiow, there exist phase-coherent 
states that extend around the physical edge of A. For 
such a state, the effect of a non-integer flux cannot be 
gauged away, lest the phase of the wavefunction become 
discontinuous at a point along the edge of A. The ener- 
gies of these edge states change as a function of $ and 
may intersect e%. Let Across label the values of flux where 
a single-particle state on the edge of A intersects e* in 
the range <I> € [0, 27r). Let n(+)($„o,J (n(_)($„o,J) be 
the number of left-edge states that cross e% with positive 
(negative) ^ at Across- Then from our definition of the 
disordered trace index, 



At^,, = lim,^o+ E [C^] (*= 



S) 



- Tr [C^] ($_ - 5) 
^ «(-)($„„3j-n(+,($c„.=) 



(21) 



The index counts the net number of chiral single-particle 
states that extend around the physical edge of A and 
whose energies intersect e*. This number is independent 
of the value of e% as long as e* remains in the mobility 
gap. The presence of these edge states at indicate a 
quantum Hall response and we identify A'^^-^^ with the 
Chern number. 

If we begin with a nontrivial Chern insulator, we 
demonstrate that this insulator undergoes a disorder- 
induced phase transition to a trivial insulator as we in- 
crease the disorder strength W. We average the trace 
index of the disordered CI over 400 random configura- 



tion of I 



with fixed ej. The disorder-quenced trace 



index is plotted as a function of W in Fig. [5] for cylin- 
ders with dimensions Nx — 10 by Ny — 5, 25 x 25 and 
40 X 40. As the dimensions of the cylinder are increased, 
we observe an increasingly sharp phase transition from 
the nontrivial insulator (with trace index 1) to the triv- 
ial insulator (with trace index 0); this occurs at around 
Wc — 5 in units t = t' = 1. 

The nonzero trace index up to the critical value of Wc 
demonstrates that the disordered CI has extended edge 
states that resist Anderson localization for moderately 
strong disorder. There remains extended bulk states that 
are spectrally connected to the edge states which allows 
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FIG. 5. Disorder-quenched trace index for the Chern insula- 
tor as a function of the disorder strength parameter W. The 
indices are plotted for cylinders of dimensions: = 10 by 
Ny = 5 (circles), 25 x 25 (crosses) and 40 x 40 (squares). We 
observe an increasingly sharp transition from the nontrivial 
insulator with (with trace index 1) to the trivial insulator 
(with trace index 0) as the dimension of the cylinder is in- 
creased. 



charge to be transported from one edge to another when 
flux is inserted - the quantum Hall effect. Hence, spectral 
flow is preserved over part of the physical spectrum. This 
conclusion is supported by studies of the level statistics 
of the physical energy spectrum and entanglement spec- 
trum of the CL^i^ 

As W increases, the range of physical energies in which 
both bulk and edge states display spectral flow is com- 
pressed and finally vanishes at Wc - this has been de- 
scribed as 'levitation and annihilation'.!^ At this point 
the edge states cease to extend around the cylinder and 
localize. The trace index vanishes because localized 
wavefunctions do not disperse with the flux and there 
are no longer any extended single-particle states that in- 
tersect in the mobility gap. 



III. QUANTUM SPIN HALL INSULATOR 

The quantum spin Hall insulator (QSHI) is a time- 
reversal invariant topological insulatop^2H56].^yjjj(,]^ j^ost 
easily modeled by combining two CI's, one for each spin 
and with opposite chiralities. This insulator is realized 
in HgTe/CdTe quantum wells and has an effective four- 
band Hamiltonian on the square lattice with periodic 
boundary conditions: 



h{k) = sin fca;ri -|- sin fc„r 
M{k) = 2 — m — cos k. 



.y.2^M{k)To 
— cos ky . 



(22) 
(23) 



where fi = cr3(8)Ti,f2 = Js(8)T2,ro = Is'S)T3. Let Is and 
lo be the identities in spin and orbital space respectively; 
cr" and T° are Pauli matrices in spin and orbital spaces re- 
spectively. The time-reversal operator is T = {i(j2®Io)K ■ 
It exhibits phases analogous to the CI i.e. it is topolog- 
ical for < m < 4 when gapped, and trivial otherwise. 
The total Chern number of the ground state vanishes 
but there is a Z2 invariant (x^^) which classifies the two 
possible equivalence classes of time-reversal symmetric 
Hamiltonians. We choose the convention that is odd 
(even) when the insulator is nontrivial (trivial). 

In our chosen geometry, we numerically plot in Fig. 
|6]-a and d the physical energy spectra of a x^^-even and 
odd insulator respectively. In the x^^.gyen trivial case, 
the conduction and valence bands are spectrally discon- 
nected; in the nontrivial case, the conduction and valence 
bands are spectrally connected by an interpolating set of 
edge states. The single-particle entanglement spectra of 
trivial (Fig. [6]-b) and nontrivial (Fig. |6]-e) insulators are 
also distinguished by a set of edge states which interpo- 
late across the entanglement gap for the nontrivial case. 
These edge states extend over the entanglement cut be- 
tween regions A and B and contribute maximally to the 
entanglement entropy. 

In principle, one may add any time-reversal symmetric 
couplings to the Hamiltonian of (22 1, including terms 
that couple the spin-up and spin-down sectors. These 
terms may change the energy spectrum and the character 
of the ground-state wavefunction. However, as long as 
both the symmetry and the energy gap are preserved, the 
topological invariant cannot change. The following 
discussion of the Z2 trace index is applicable to the most 
general form of a time-reversal invariant Hamiltonian. 



A. QSHI: The trace index and its equivalence to 
the Z2 invariant 

In this section we define the trace index for a Z2 insu- 
lator and show its equivalence to the Z2 invariant. The 
differences from the U{1) trace index are that the Z2 in- 
dex is defined over half the Brillouin zone and only the 
parity of the index is physically relevant. 

The correlation matrix G(ky) as defined in (j3| satis- 
fies TG{ky)T~^ = G{—ky) if the subspace of occupied 
bands at ±fcy in the physical spectrum is closed under 
time-reversal. This condition means that if any single- 
particle state \ip) is occupied, so is its time-reversed part- 
ner T\'ijj). We prove this claim in Appendix |V] We as- 
sume this is true for all ky and will later generalize our 
discussion. The reduced correlation matrix G^{ky) as 
defined in (j4]) is the matrix G{ky) with spatial indices 
restricted to region A. Since the time-reversal operation 
T — K ZCT2 ® lo is diagonal in spatial indices, a spatial 
entanglement cut cannot break time-reversal symmetry 
and TG^{ky)T-^ = G^{-ky) as well. This symmetry al- 
lows us to extract all relevant information from half the 
Brillouin zone {ky G [0,7r]). 
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FIG. 6. Top: Figures [6]-a,b,c show, respectively, the physical energy spectrum, the single-particle entanglement spectrum (ES), 
and the trace of G^(ky) for a trivial insulator (x^^ even) on an open-end cylinder. In Fig. [Gl-a and b, bulk bands on opposite 
sides of the gap are spectrally disconnected from each other. In Fig. [6|c, the trace of G is a continuous function of ky so 
the trace index Az.2 = 0. Top: Figures |6]-d,e,f show, respectively, the physical energy spectrum, the single-particle ES, and the 
trace of G^{ky) for a nontrivial insulator with odd ■ In Fig. [sj-d and e, a set of edge states interpolate between bulk bands 
on opposite sides of the gap - this spectral flow distinguishes trivial and nontrivial Z2 insulators. In Fig. [oj-f, the trace of 
jumps discontinuously by 1 in half the Brillouin zone ky £ [0,7r] - the trace index Az2 = +1- 



The following discussion mirrors that in Section II A 



Let Kl^l^^ be the values of ky in half the Brillouin zone 
where states that are localized on the edge of A intersect 
the Fermi level. Our first example is a x^^-odd insula- 
tor with one such intersection at K^l^^J — As illus- 
trated in Fig. [7]-a, its physical energy spectrum exhibits 
spectral flow - the conduction and valence bulk bands 
are spectrally connected across the gap by an interpo- 
lating set of edge states. We choose our Fermi energy 
to lie just above the Kramers' doublet of left-edge states 
at ky — 0, so each occupied single-particle state has a 
Kramers' partner which is also occupied. Following the 
discussion in the previous paragraph, we conclude that 
TG'^{ky)T-^ = G^{~ky) for all ky. We plot the time- 
reversal symmetric spectra of in Fig. [7]-c. 

We define the Z2 trace index (Az^) for clean, free in- 
sulators as the sum of the trace discontinuities over all 
^cioaJ ill the thermodynamic limit: 



Az2 = lim V Tr [G^] (K^:^^ + ^) 



-Tr {G^\{K'-^J2 - w) 



(24) 



As illustrated in Fig. [7]-e, the trace of G^ decreases dis- 
continuously by 1 between ky = and in our example, 
so = — 1 and has the same parity as x^^- 



More generally, we claim that Az^ is modulo 2 the to- 
tal number of left-edge states that intersect the Fermi 
level in half the Brillouin zone. Following the discus- 
sion in Section |II A| we apply Theorem |II.1| and obtain 



Az2 — 



where 



is the number of 



left-edge states that intersect the Fermi energy with pos- 
itive (negative) . The claim is proven by the identity 

n(_) — n(+) = -|- n^^f mod 2. 

If n(_i -t- is odd, the physical energy spectrum 
is characterized by spectral flow, which cannot be re- 
moved by any change in the Hamiltonian that simultane- 
ously preserves the gap and time-reversal symmetry - this 
phase displays topologically protected edge states and 
has odd x^2_ The trivial phase corresponds to both x^^ 



and ' 



being even, in which case there is no spec- 



tral flow. Since Az2 = n^. 



mod 2, we identify the 



Z2 trace index with the Z2 invariant - Az2 = X^^ mod 2 
- with the caveat that the subspace of occupied bands 
in the many-body ground state is closed under time- 
reversal. We numerically plot the trace of G^ for the 
trivial and nontrivial insulators in Fig. |6]-c and f respec- 
tively. The Fermi level of the nontrivial ground-state is 
depicted in Fig. [6]-d as a dashed line; the subspace of oc- 
cupied bands is closed under time-reversal. In the trivial 
case, Tr G"* is a continuous function of ky while in the 
nontrivial case it has a discontinuity of -1-1 in half the 
Brillouin zone, so we confirm that the parities of 
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FIG. 7. Left side: Figures [7]-a,c,e show scliematically the physical energy spectrum, the single-particle entanglement spectrum 
(ES) and the trace of the single-particle ES respectively, as functions of ky, for a many-body ground state of a nontrivial QSHI. 
The subspace of occupied bands in the physical spectrum is closed under time-reversal. Right-side: Figures [7}b,d,f label the 
new functions when an integer charge-flux is inserted through the symmetry axis of the cylinder. For convenience, bulk bands 
away from the gaps are portrayed as dispersionless. Eigenvalues that are colored red, green and blue correspond to eigenstates 
that extend over the left edge, the entanglement cut and the right edge of the cylinder respectively. Kramers' doublets at 
ky = and vr are labelled by thick dashes in the physical energy spectra; cccupied single-particle states of a many-body Slater 
determinant are labelled by black circles. 



and are identical. 

We now address the exceptional case of a ground state 
for which the subspace of occupied bands at the sym- 
metric momenta fc^"™ = 0, tt is not closed under time- 
reversal. If at ky''" the Kramers' doublet that extends 
around the edge of A is singly-occupied, the single- 
particle entanglement spectrum (ES) is no longer time- 
reversal symmetric, i.e. TG^{k^''"')T^^ ^ 
This is apparent because a singly-occupied doublet con- 
tributes one zero eigenvalue and one unit eigenvalue to 
the single-particle ES at fc^"", in accordance with Theo- 
rem |II.1[ and neither of these eigenvalues has a degener- 
ate Kramers' partner. 

As an example, we consider an insulator with an odd 
number of particles, which implies that the ground state 
is two-fold degenerate. Our previous example of a many- 
body ground state in Eig. [7] has an even number of par- 
ticles, so we remove a particle \^) from the two-fold de- 



generate subspace at energy Ep and momentum ky = 0. 
This state {ip) may be any linear combination of the ba- 
sis states that span the Kramer's subspace at ky — 0. 
Since = — /, the time-reversed state T lip) is always 
orthogonal to and hence lies within the occupied sub- 
space - we have thus shown that the occupied subspace 
is not closed under the action of time-reversal. As illus- 
trated in Eig. |8]-a, the resultant Slater determinant has 
a singly-occupied left-edge Kramers' doublet at ky — 0. 
The occupied left-edge state in the physical spectrum 
translates to an eigenvalue of 1 at /cy = in the single- 
particle ES of Eig. [8]-c; its unoccupied Kramer's partner 
in the physical spectrum translates to an eigenvalue of 
in Eig. |8]-c. There are now two isolated single-particle 
entanglement states, separated across the entanglement 
gap, at symmetric momentum ky — 0; neither of these 
two states has a Kramer's partner. Since there is no dis- 
continuous change in the occupation numbers of left-edge 
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FIG. 8. Left side: Figures [8]-a,c,e show schematically the physical energy spectrum, the single-particle entanglement spectrum 
(ES) and the entanglement energy spectrum respectively, as functions of ky, for a many-body ground state of a nontrivial QSHI. 
The subspace of occupied bands in the physical spectrum is not closed under time-reversal. Right-side: Figures [8]-b,d,f label 
the new functions after tt charge-flux and n F-flux are inserted in succession through the symmetry axis of the cylinder. For 
convenience, bulk bands away from the gaps are portrayed as dispersionless. Eigenvalues that are colored red, green and blue 
correspond to eigenstates that extend over the left edge, the entanglement cut and the right edge of the cylinder respectively. 
Kramers' doublets at fc^ = and n are labelled by thick dashes in the physical and entanglement energy spectra; occupied 
single-particle states of a many-body Slater determinant are labelled by black circles. 



states, the trace of G"* is now continuous over the half 
Brillouin zone and the Z2 trace index is zero. We have 
thus shown that the parity of is opposite to that of 
for an exceptional case. 
We remedy this index by introducing a quantity Q 
which we define as the total number of singly-degenerate 
entanglement states above the entanglement gap at both 
ky^"' . In Appendix VI we remark briefly on how to rig- 
orously define Q for even and odd Ny. Adding Q to 



gives the correct parity of for any value of the Fermi 
energy: 



Az2+<d^ X^^ mod 2. 



(25) 



For our exceptional example, Fig. [8]-c reveals a singly- 
degenerate entanglement eigenstate with unit eigenvalue 
above the entanglement gap, so 8 = 1 and Az., + B gives 
the correct parity of ■ 
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B. QSHI: Flux threading and entanglement 



In this section wc investigate the effects of threading flux through the symmetry axis of a QSHI in an open-end 
cyhnder geometry. Unhke the CI, the QSHI is constrained by time-reversal symmetry to have zero Chern number 
and it does not pump any net charge in response to charge flux - flux that induces a spin- and orbital-independent 
response. For the QSHI, we also consider flux that differentiates electrons with different spin and/or orbitals. The 
flux threading may be implemented by modifying the hoppings between sites with y = 1 and y = Ny so that they 
carry a Peierls's matrix g*27rr*/$o -v^ith T a Hermitian matrix in spin-orbital space - we call this matrix or F-flux. The 
tight-binding model in ([T]) now has the nearest-neighbor form 



xy,ap 



— ^ '^xyaf^(x,y){x,y)^xyP ~^ [^l+lya^{x+l,y){x,y)'^xyP ~^ 4cy+laf^'(x,y+l){x,y)[' 



<I> - 



Isp^xy/S 



+ h.c. (26) 



with the matrix h"^^ obtained from Fourier transforming the QSHI model (22) in i and y. a, (5 are spin/orbital indices 
and i, j, {x,y) are lattice coordinates. If the matrices ^(ia;)(Ar^£c) 



and e*^'^ do not commute, the Hamiltonian ( 26 ) 



cannot be reduced to a flux-free Hamiltonian by a different choice of gauge. In this case, the energy spectrum is 
changed by inserting matrix flux and translational symmetry along y is lost. 

I 



In Section |III B 1| we investigate the effect of thread- 
ing charge flux {T = Ig (E) lo) on the trace index of the 
single-particle entanglement spectrum. We show that a 
certain choice of F induces charge transport in Section 
IIIIB 21 and demonstrate another method to extract the 
Z2 invariant. In Section IIIB3 we identify the signa- 
ture of charge transport in the many-body entanglement 
spectrum as a flow in particle-number space. 



1. Charge flux and time-reversal polarization: effect on 
entanglement 



Because the QSHI Hamiltonian has global U{1) sym- 
metry and T = Is ® lo commutes with ^(i2;)(Arj,x)i the 
charge-flux-inserted Hamiltonian is gauge-equivalent to 
a flux-free Hamiltonian. The physical energy spectrum 
is thus invariant and ky is still a good quantum number. 
The effect of a unit charge flux is to transform single- 
particle states with ky to adjacent states at ky + j^. 
Unlike the CI, there is no net transfer of charge be- 
tween edges because every inter-edge current has a time- 
reversed partner. Nevertheless, the occupation numbers 
of the single-particle states change due to a quantum spin 
Hall response. We illustrate this process for a x^^-odd in- 
sulator with a many-body ground state depicted in Fig. 
[7]-a; after inserting a unit flux, the corresponding flux- 
evolved state is shown in Fig. [7[b. The resultant wave- 
function is [7[]^7"[7ij]^7ij^ IV'gs) where 7^^ and 7^ are 
normal mode operators of the lowest energy states on the 
right edge and /, // denote time-reversed partner modes. 

As with the Chern insulator, the positions of the dis- 
continuities in the trace are translated by K'^l^^J — > 



2tt 



but the magnitudes of the discontinuities 



are invariant. Since the Z2 trace index is defined 
over half the Brillouin zone {ky £ [0, tt]), it is possible to 
(i) translate a discontinuity at KlU^^ ~ ~ '-'^^ 
the half-BZ or (ii) translate a discontinuity at — into 



the half-BZ. For our example, the discontinuity of -f 1 at 
— in Fig. j7|e (before) is translated into the half-BZ 
at -t-j^ in Fig. [7|f (after). Before, Az„ = -1; after, the 
discontinuities at 



and cancel so Az^ — 0. 



We have shown that the index Az2 is not invariant un- 
der flux insertion. However, (i) is always accompanied by 
a change in parity of the occupancy of the ky — tt left- 
edge Kramers' doublet (in the physical spectrum) and 
similarly (ii) changes the parity at ky = 0. The resul- 
tant change in parity of O in the entanglement spectrum 
compensates the change in parity of Az^ so that the net 
parity is invariant. For our example, we may compare 
the spectra of in Fig. [7]-c (before flux) and d (after). 
The difference is that the flux-evolved spectra of G'^ has 
an additional non-degenerate eigenvalue of 1 at fcj^ = 0, 
thus flux-insertion has changed the parity of Q. 

In Section HI A we identified Az2 + © with the Z2 
invariant of the ground state; in this Section we show 
that this equality may be extended to any flux-inserted 
state. This implies that we may separately extract 
from both the ground-state and its flux-inserted state 
through Azn + 0- 



2. Matrix flux and charge transport: effect on entanglement 

In a generalization of the Hall effect for Z2 insulators, 
Qi and Zhang proposed that if F is odd under time- 
reversal (TFT-i = -F) and satisfles e"^ = in- 
serting the matrix flux e»27r*/*or i]^(jm;es charge trans- 
port between edges. ''-^ Henceforth, we call flux that sat- 
isfies these conditions as F-flux. The first condition en- 
sures that the Peierls's matrix satisfies T e^^^r^/^o j^-i _ 
gi2irr*/$o go the Peierls-modified Hamiltonian is time- 
reversal symmetric for all values of The second condi- 
tion implies that for special values of the flux ($ = n$o /2 
with n an integer), the matrix g«2:iTr$/'i>o proportional 
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to the identity and commutes with h 



iix)iNyxy, TT (27r) 
F-fiux has the same effect as tt (27r) charge flux on the 
physical spectrum. 

A choice of F that satisfies both conditions is fiiai (E) 
lo where n is a unit vector. This choice is caUed spin 
flux because it picks a preferred direction (fi) in spin 
space. Another vahd choice, T = Is®T2, picks a preferred 
direction in orbital space. The authors in Ref. ,79, showed 
that different choices of F cannot change the parity of the 
charge pumped during an adiabatic insertion of tt of F- 
flux; this parity is precisely the Z2 invariant ■ 

We now present a way to extract from two many- 
body states differing by tt F-flux through the entan- 
glement spectrum; this method applies to disordered 
and/or interacting Z2 insulators as well. In Section IIB 1 
we showed that the full trace of the reduced correla- 
tion matrix is the number of particles in region A, i.e. 
Tr = {^NaY^)- Let and G^^ respectively de- 
note the reduced correlation matrices before and after 
insertion of tt F-flux. The difference in the trace of these 
two matrices is the change in the number of particles 
in region A, which is also the charged pumped when tt 
F-flux is threaded: 



Tr Gtr - Tr = x''" mod 2. 



(27) 



In Fig. [9} a we thread spin flux (F = (Ti (g) /q) sepa- 
rately through a trivial and nontrivial disordered QSHI 
and plot the trace of G"^ as a function of flux $ for each 
case. The result is that Tr G^^ — TrG^ — l for the non- 
trivial insulator with odd and equals for the trivial 
insulator with even x^^- If instead we thread orbital 
flux (F = Js (g) T2), we show in Fig. [9]-b that no charge is 
pumped for the trivial insulator and Tr G^p— Tr G^ = — 1 
for the nontrivial insulator. In either case, the parities of 
the pumped charges do not change with a different choice 
of F. 



3. Spectral flow and the many-body entanglement spectrum 

In this section we investigate how spectral flow man- 
ifests itself in the many-body entanglement spectrum 
(ES) of a Z2 insulator with a cylinder geometry. Since 
e^'^^ = —I, inserting tt F-flux has the same effect on the 
physical spectrum as inserting tt charge flux. To make a 
closed loop in parameter space we thread in succession tt 
F-flux and tt charge flux (ttTttI). Since the Z2 insulator 
is time-reversal symmetric, charge is pumped only in the 
flrst half of the loop when we thread F-flux. The parity 
of the net charge that is pumped in this closed loop is 
thus the Z2 invariant. Let us deflne the amount of charge 
pumped from the right to left edge in this adiabatic loop 
by Qr- The absolute value of Qr depends on the choice 
of F, however we have shown that Qr — X^^ mod 2. 

The following discussion is analogous to the discussion 
of the CI in Sec. 



IIB2 



Let 1*0) and \'^^r^i) be the 
many-body wavefunctions before and after the insertion 
respectively. The physical energy spectrum is invariant 
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FIG. 9. Trace of vs flux $ for a trivial x^^-even insula- 
tor (red) and a nontrivial x^^-odd insulator (blue). Flux is 
threaded by inserting a matrix e'2irr#/*o nearest-neighbor 
hoppings between y = 10 to y — 1 on a 10 x 10 lattice. 
Our Hamiltonian is Eq. ( 22 1 with (i) an additional inversion- 
symmetry-breaking term 0.1 sinfca;(Ti ® /□ + 0.1 smkya2 ® T3 
that couples spin-up and spin-down sectors and (ii) a weak 
random on-site potential so that the insulator is disordered; 



X 



is odd (even) when m = 1 (— 1). In Fig. [9]-a (b) we choose 
V = G\® lo {la ® T2). For both choices of F, the trace of 
smoothly changes by 1 as tt F-flux is threaded through the 
nontrivial insulator. In the trivial case, the traces at tt and 
flux are equal. 



under TrFTrJ flux but the number of occupied edge states 
in region A increases by Qr- As an example we thread 
ttTttI flux into the many-body ground state of a x^^- 
odd insulator as depicted in Fig. |8]-a. We may pick F 
such that Qr = 1; the resultant wavefunction is shown 
in Fig. |8]-b to have an additional left-edge particle at 
ky = j^. We compare the entanglement energy spec- 
tra before (Fig. [sj-e) and after ([sj-f) flux threading; the 
difference is that the flux-excited insulator has an addi- 
tional eigenvalue ^ = — cxd at fcj, = This means that 
if we began with an entanglement ground state with 
particles under a Fermi level ^j?, the flux-evolved entan- 
glement ground state now has + 1 particles under the 
same Fermi level In Sec. |IIB 2 we argued that the 
extra entanglement energy at ^ = — 00 cannot change the 
entanglement entropy, so the sector of the many-body 
ES before threading is identical to the + 1 sector after 
threading. 

More generally, we apply Theorem |II.1| and the proof 
in Appendix |IV| to derive a flow between particle number 
sectors: 



1*0) ( -) I'S^r^i) 
Na J Na+Qy 



(28) 



This is the signature of spectral flow in the many-body ES 
of a Z2 insulator. For interacting insulators that preserve 
the global U{\) symmetry, this relation is valid as well. 
We note that the parity of the number sector is changed 
if the QSHI is nontrivial and unchanged if trivial. 

Let us deflne the even (odd) particle-number sector as 
the set of non-zero eigenvalues of Pa obtained by occu- 
pying the single-particle entanglement eigenstates Ci of 
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Eq. ([8]) in all possible combinations such that the total 
occupation number is even (odd): 



even (odd) 



Na 



Na even (odd) 



(29) 



As a corollary to (28), we deduce that inserting nTnl 



flux through a nontrivial Z2 insulator results in the even 
number sector flowing into the odd number sector (and 
vice versa), so 



1*0) 



X^2odd 



even (odd) 



A . (30) 

odd (even) 



This is distinguished from the trivial Z2 insulator, for 
which the parity of the number sector is invariant under 
the same flux insertion. 



IV. CONCLUSION 

The fundamental property that distinguishes free- 
fermion topological insulators which produce a quantum 
Hall response is spectral flow, which is manifest in both 
the physical energy spectrum and the single-particle en- 
tanglement spectrum (ES) for a spatial entanglement cut. 
By a quantum Hall response, we refer to both time- 
reversal breaking Chern insulators and time-reversal in- 
variant Z2 insulators. For an open-end Laughlin cylinder 
geometry, spectral flow is identified by the presence of a 
set of edge states that interpolate across the energy gap 
in both physical and entanglement spectra, such that the 
conduction and valence bands are spectrally connected. 

In principle, there are many ways to partition the 
Hilbcrt space in an entanglement cut; we would like a 



cut that explicitly demonstrates the connection between 
entanglement and quantum Hall physics. In this article 
we have found that a spatial entanglement cut that sep- 
arates the two edges of the cylinder is ideally suited to 
probing the quantum Hall response of a topological in- 
sulator. During the adiabatic insertion of an appropriate 
flux, the quantum Hall current flows directly from one 
edge of the cylinder to the opposite edge, or equivalently 
from one subsystem of the entanglement cut to the other. 
Since the quantities of pumped charge for both Chern 
and Z2 insulators are directly related to their respective 
topological invariants, it must be that these topological 
invariants are encoded in the ES. 

In this work we have shown that the signatures of 
spectral flow in the ES lie in (i) a non-vanishing 'trace 
index' of the single-particle ES (ii) the adiabatic evolu- 
tion of the full trace of the single-particle ES when flux is 
threaded (iii) a global translation in the particle-number 
sectors of the many-body ES when flux is threaded. 
Any of the three signatures provides a means to extract 
the desired topological invariants. While all three are 
generalizable to disordered insulators, only (ii) and (iii) 
are applicable to interacting insulators. 
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Appendix A: Derivation of Eq. [T] 



The one-body reduced correlation matrix as defined in Q is 



ql 



11 1-1- £-«,(*:„) ^kyqXkyl 



(Al) 
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Let \4'ky) = \^iky ■ ■ ■ n-QAky) ^Sibel a Slater determinant in the momentum sector ky with the corresponding occupation 
numbers {niky}. We need the identity 



Tr.fc., 



i=l 
1 



E ('^'^■.(1)1 n 



ihyi^)) (0fey(2)| xLqXkyl I'l^kyil)) 



J2 e-«-('=«)"-«(2) (^fc„(l)| ^fe^(2)) (0fe„(2)| xI„,Xfc„; 



1 

n E ^ 

i=l nifc„=0 



{niky 



Qa 1 

i=l "ifc =0 



2— 



The second equahty foUows because e 



is a counting operator that is diagonal in the basis 1 4>ky ) • 



(A2) 



Substitution of this identity in ( Al ) reveals 



Since it*^^^ is confined in A, uf' = if j G B. This means 
that 



(3) 



Recall that the entanglement Hamiltonian matrix is 

[i?.„,]^^.(fc,) = J2 <^^v) ^T^^v) ^^^^y) (4) 

9 

Hence Eq. 

II. PROOF OF THEOREM ITlT] 



Vie A, 
Vie A, 



(L)* 



(3) 



Hence we have shown that u^-^^* is an eigenvector of G'^ 
with zero eigenvalue. By definition, 



(^A _ ^A I p{L)A* 

ij ij ij 



Gt,+v^^*v!t^ (4) 



Proof From Eq. (3|, we know that the complex con- 
jugate of the correlation matrix G is a sum of projec- 
tors of single-particle eigenstates of the Hamiltonian H p(L)A* ^j^g projector of the eigenstate u^^^* which we 
which are occupied in the many-body wavefunction. Let have shown to be a zero-eigenvalue eigenvector of G^. In 
us define G as the correlation matrix of a many-body other words, w'^^* is orthogonal to all the eigenvectors 



wavefunction that includes the single-particle state u 



of G'*. Since G is Hermitian, so is G^; this implies that 



Since w'^^ is an eigenstate of H we may form a projector G'* is orthogonally diagonalizable. If G^ is an n x n ma- 
matrix p/f ^ = ui^^'u-^^* and define G as the difference trix, we can find n linearly independent and orthonormal 

eigenvectors v that diagonalizes G'*. Let A*^*^ denote the 
eigenvalues of G^ with corresponding eigenvector v'^^\ 
then 



between G and F^^'* 



G = G - P^-^)* 



(1) 



G is the correlation matrix of a many-body wavefunction 
that excludes the single-particle state w'^). We define G^ 
and p(^)"^ respectively as the matrices G and p(^' with 
indices restricted to A. Since G* is a sum of projectors 
of states orthogonal to u^^\ 



jeA(giB 



(2) 



G' 



(5) 



9=2 



with w'-'^^ orthogonal to the matrix J2i=2 A'-'-'f In 
fact, we have already found such a form for G"* in Eq. (|4[ 
and we need only identify the eigenvector w*^^^ as 
and the eigenvalue A'^^ as unity. | 
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III. PROOF OF INVARIANCE OF THE U{1) 
TRACE INDEX UNDER INTEGER FLUX 
INSERTION 

For a generic CI with Chern number Ci, the left-edge 



states intersect the Fermi energy ri+ times with 



de 



> 



and 



times with 



dk^j 



< 0, such that 



Ci . There are n+ values of K^^_,^ where the trace of 
discontinuously decreases by 1 and n_ values of K~^^^ 
where the trace of G"* discontinuously increases by 1. 

for each there is at 



II.l 



According to Theorem 
least one occupied left-edge state at K^^, 

is a unit-eigenvalue eigenstate of G^{K^^^^ ~ ^) ''^^'^ 

least one unoccupied left-edge state at K'^ , + which 

cross JVy 

is a zero-eigenvalue eigenstate of G"^(K'^_ , + ^) - each 
of these crossings contributes —1 to the trace index. For 
all Kz„.. there is at least one unoccupied left-edge state 

— which is a zero-eigenvalue eigenstate of 

— and at least one occupied left-edge state 
+ which is a unit-eigenvalue eigenstate of 

+ - each such crossing contributes +1 to 
the trace index. 

Upon threading a unit flux, all states transform from 
momentum ky to ky + j^. We assume that the flux 
is small enough that it does not change the occupation 
numbers of bulk extended states, to keep the phase insu- 
lating. For each K^^^^, the left-edge state at K^^^^ + 

is now occupied and is a unit-eigenvalue eigenstate of 
G^iK'^ , -I- T^). In the same left-edge mode, the left- 

edge state at i^^toss + |r ^^i^^ unoccupied and is a 
zero-eigenvalue eigenstate of G^{K'^^^^_^ + If")- ^'-'^ each 
the trace of G'^ is now continuous between k,. 



at K-^^, 
at KZ,,, 



K+ 



N. 



there is now a discontinuity 



of —1 in the trace of G^ between ky 



K+ 



N„ 



and 



K 



+ 



37r 



there is 
between 



Analogously for the negative- crossings, 
now a discontinuity of -1-1 in the trace of G"* 
^y = ^c-ss + ]f^andi^-_ + ^. 

We have shown that the net result of integer-flux in- 
sertion is that all trace discontinuities are translated by 
1^. The U(X) trace index is the sum of all these dis- 
continuities and is therefore invariant under integer flux 
insertion. 



the number of every other eigenvalue, which we label by 
Cf^". We define ri'i''^^ as the occupation numbers of the 
single-particle particle eigenstates with eigenvalues Cl'^^ ■ 
By definition, d = tiq + 't-i + 't-i/2- 

All eigenvalues A in the many-body entanglement spec- 
trum are proportional to the factor Hie unoccupied (1 ~ Ci) 
- any nonzero eigenvalue A must arise from occupying all 
ni single-particle states with eigenvalue 1, so the mini- 
mum particle-number sector is ni. 

A is also proportional to the factor Ilieoccupicd Ci- This 
means we cannot have a particle-number sector larger 
than 7ii/2 -l- TLi because that would involve occupying at 
least one zero-eigenvalue single-particle state. 

For a number sector ni + m which falls in the range 
bounded by ni and ni +rii/2, we have from Eq. (|9| that 



= n c n (1-0) E"'="i 

ni+m I i—i j — i 



(1) 



with the i (j) labelling occupied (unoccupied) entan- 
glement states. For A to be nonzero, all ni unit- 
eigenvalue single-particle states are occupied and all no 
zero-eigenvalue states are unoccupied. This results in 



ni+m ni m rr, 

i=l j = l fe=l fe=l 



/2) 



(2) 



for occupied states and 



-n\—nQ — m 



d—ni—m no 

n (1-0)- n (1-0) n a-cr 

i = l k = l 1 = 1 



/2)^ 



n 1 / 9 — m 



n (i-cr') 



(3) 



1=1 



for unoccupied states. We substitute Eq. ^ and ([3| into 
([!])- up to exponential finite-size corrections, the set of 
nonzero eigenvalues in the particle- number sector rii+m 
is 



"1/2-™ 



rti +m 



in ' "1/ z — 

ncr' n (i-cr 
i=i j=i 

E(l/2) I 
nl —my. 



(4) 



IV. PROOF OF EQ. (18) AND (281 



According to Theorem |II.1[ occupied left-edge states 
are eigenstates of the reduced correlation matrix G'* with 
eigenvalue 1 with finite-size exponential corrections while 
unoccupied left-edge states have eigenvalue 0. We define 
no (ni) as the number of eigenvalues that flow exponen- 
tially to 0(1) in the thermodynamic limit, and rii/2 as 



Once again, i (j) labels occupied (unoccupied) entan- 
glement states. There are (",\{^) combinations of single- 
particle occupation numbers that give nonzero values of 
A in the ni + m sector. 

We observe that the set of nonzero eigenvalues in the 
ni + m sector is independent of ni. This result is applica- 
ble to both CI and QSHI. By inducing charge transport 
in either CI or QSHI through fiux threading, the wave- 
function only changes in the occupation numbers of edge 
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states, hence by Thm. II. 1 only ni is affected, while Ci^''^' 
are invariant. 

Suppose we began with a = ni+m number sector. 
If Q charges are transported from the right to the left 
edge due to flux threading, ni increases by Q but the set 
of nonzero eigenvalues is invariant. The set of many-body 
entanglement eigenvalues in the Na sector before flux 
threading is identical to the set of eigenvalues in the + 
Q sector after flux threading. This statement is expressed 
mathematically in Eq. (18) and (28) in the contexts of 



the CI and QSHI respectively. The statement is true for 
all Na so there is a global translation in particle-number 
space. 



V. PROOF THAT TG{ky)T-^ = G{-ky) IF THE 
OCCUPIED SUBSPACE AT ±ky IS CLOSED 
UNDER TIME-REVERSAL 

We write the time-reveral operation as T — V K with 
V a unitary operator and K complex conjugation. T 
acts on spin-half single-particle states as = —I. We 
define the time-reversal sewing matrix in the basis of oc- 
cupied bands Fnm{ky) — {u'^{—ky)\ T \u^{ky)) ; 71,111 — 
1 . . . N^^^ and another sewing matrix in the basis of 
afl bands Hnmiky) = {u"{—ky)\T\u"^{ky)); n,m — 
1 . . . N^ii . F is unitary if the subspace of occupied bands 
is closed under time-reversal, i.e. H decouples into two 
block diagonals for the occupied and unoccupied bands. 
Proof: Using the completeness property of the Bloch 
eigenvectors and the summation convention, we show 
that H is unitary: 



Hjiiky)Hjm{ky) 

- <{-ky) Uliky) Ui*i-ky) Vj^*{ky) = 5„„ 



If H decouples into two block diagonals for occupied and 
unoccupied bands, the unoccupied block F must itself be 
unitary. If F is unitary, we derive two identities from the 
definition of F: 

V„U Uili-ky) = Fij{ky) Uj{ky) 

Fuiky) Ul,{-ky) = Vran u\* {ky) . (2) 

From these two identities, V = ia2 <E) lo — V* , and the 
definition of G{ky) in Eq. (|3]), 

iTG{ky)u = Km^<:(fc,)"f(fc,) 

= <iky) U]*{ky) K = F:^{ky) <(-fc,) U\{ky) K 

= Ki^ky) uU-ky) K = {G{-~ky) TU 
TG{ky)T-^ =G{-ky) as advertised. (3) 



VI. DEFINING e FOR EVEN AND ODD Ny 

We assume periodic boundary conditions unless flux 
is inserted. If Ny is even, ky — Q and ky — n are al- 
lowed momenta for physical states, and we define — 



^ 0;ky = 0) 



)(<& = 0;ky = tt) as the sum 



of singly-degenerate entanglement eigenstates above the 
entanglement gap at = and tt. 

If Ny is odd, fcy = is an allowed momentum for 
physical states but ky = tt is not. However, upon 
adiabatic insertion of tt charge flux we now have an- 
tiperiodic boundary conditions, and the single-particle 



states at fc„ 



flow to TT. 



Now fcj^ = TT is al- 



(1) 



lowed but ky = is not. Hence we define for odd Ny, 
e = n(,)($ = 0;ky = 0)-t-n(,,($ = ^;fcy = tt) where the 
first term is evaluated with periodic boundary conditions 
and the second with antiperiodic boundary conditions. 
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